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ABSTRACT 

We evaluate the ability of the BATSE spectroscopy detectors to detect the 
absorption lines observed by Ginga in gamma ray burst spectra. We find that 
BATSE can detect the 20.6 keV line in the SI segment of GB870303 with a 
detection probability of ~ 1/4 in nearly normally incident bursts, with the 
probability dropping off to nearly at a burst angle of 50°; the lines at 19.4 
and 38.8 keV lines in GB880205 have a high detection probability in BATSE 
up to burst angles of 75°. In addition, we calculate detection probabilities for 
these two line types as a function of signal-to-noise ratio and burst angle for 
use in detailed comparisons between BATSE and Ginga. Finally, we consider 
the probability averaged over the sky of detecting a line feature with the actual 
array of BATSE detectors on GGRO. 

Subject headings: gamma rays: bursts 



1. INTRODUCTION 

What are the line detecting capabilities of the Burst and Transient Source Experiment 
(BATSE) on board the Gompton Gamma Ray Observatory (GGR0)1 As described by the 
first paper in this series (Palmer et al. 1994b), currently no absorption features have been 
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detected in the spectra observed by BATSE's spectroscopy detectors (SDs), while two sets 
of hnes were discovered in the bursts observed by the Ginga gamma ray burst detector. 
In the second paper of this series (Band et aL 1994) we presented both Bayesian and 
"frequentist" methodologies for comparing the consistency of the Ginga and BATSE results. 
These comparisons require quantitative detection probabilities for lines in each spectrum 
accumulated by the two detectors; in the present paper we calculate these probabilities 
for the BATSE detectors. In addition, if hnes are detected in the BATSE bursts, we 
will then estimate the rate at which lines occur, which also requires the line detectability 
probabilities. 

The existence of absorption lines is one of the crucial remaining questions which 
BATSE can address in the study of GRBs. One of the key pieces of evidence supporting the 
Galactic neutron star origin of GRBs was the series of reported absorption features in the 
15-75 keV band (Konus— Mazets et al. 1981; HEAO i— Hueter 1987; Gm^a— Murakami et 
al. 1988) interpreted as cyclotron absorption in a teragauss magnetic field (e.g., Wang et al. 
1989); the only known astrophysical body with such strong fields is a neutron star. BATSE 
found that the burst distribution is isotropic yet radially non-uniform (Meegan ct al. 1992), 
overturning the Galactic disk neutron star paradigm; whether BATSE also finds absorption 
lines is the subject of this series of papers. If these lines are discovered they will be an 
important probe of the burst environment; in particular, it will be difficult (but probably 
not impossible) for cosmological models to explain absorption features. However, in the 
absence of any conclusive BATSE detections (Teegarden et al. 1993; Band et al. 1993a; 
Palmer et al. 1993, 1994a,b) we have to consider the consistency between the BATSE and 
Ginga results (Band ct al. 1994), the calculation of which requires the number of bursts 
observed by each instrument in which lines would be detectable. 

To properly study the rate at which lines occur in the bursts observed by past and 
present missions, we need both a quantitative assessment of the reality of the claimed lines 
(as described below) and a list of the spectral parameters (e.g., the signal-to- noise ratio) 
which determine line detectability. Unfortunately, fit parameters were not provided for the 
line features reported in ~ 20% of the Konus bursts (Mazets et al. 1981); in addition, the 
searched spectra have not been described in quantitative detail. The two line detections 
reported among the 21 HEAO 1 bursts (Hueter 1987) are only of moderate significance (by 
the methods described below) and few details have been provided concerning the spectra 
searched in this data set. However, the burst data necessary to evaluate both the fine 
candidates and line detectabihty in the bursts observed by BATSE and Ginga are available. 
Consequently our quantitative comparison is between these detectors. 

The absorption features are undoubtedly characterized by distributions of line 
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parameters. However, there have been too few definitive fine detections to determine 
these distributions. Therefore we use the two Ginga detections as exemplars of the 
hnes which might be present. While four sets of lines were reported, the sets in the S2 
segment of GB870303 (Murakami et al. 1988) and in GB890929 (Yoshida ct al. 1991) 
are not sufficiently significant (by the detection criteria discussed below) to be considered 
detections, while the lines in the SI segment of GB870303 (Graziani et al. 1992) and in 
GB880205 (Murakami et al. 1988) are regarded as real. Therefore we study the detectabihty 
of the single line in the SI segment of GB870303 (20.6 keV) and the harmonically spaced 
hnes in GB880205 (19.4 and 38.8 keV) in the BATSE detectors as functions of the relevant 
parameters. 

In this paper we explore the detectabihty of absorption lines as a function of various 
parameters by the BATSE SDs, both to demonstrate that BATSE can indeed detect lines 
similar to those observed by Ginga, and to derive detection probabilities which will be 
useful for future studies. After describing the detector (§2), we outline the methodology 
used for this study (§3). The presentation of our results follows (§4). Finally, we analyze 
the dependence of the detectabihty on various instrumental parameters, and consider the 
average detectabihty of Gm^'a-like lines in the BATSE detectors (§5). The implications of 
these detection probabilities for BATSE- Gm^'a consistency will be discussed later in this 
series of papers. 



2. THE DETECTORS 

Because more complete descriptions of the BATSE SDs (Fishman et al. 1989, 1995; 
Band et al. 1992) are available elsewhere, here we provide a brief overview of the detectors, 
emphasizing aspects which are important for understanding their line detection capabilities. 

Each of BATSE's eight modules contains an SD, and therefore every burst is observed 
by a number of these detectors; however, not all the SDs from which burst data are available 
may have been configured to provide useful spectra (e.g., as a result of the discriminator 

and gain settings). The SDs are simple scintillation detectors built around a 5" diameter 
by 3" thick cylindrical Nal(Tl) crystal with no active shielding. The detectors are sensitive 
to gamma rays entering the side of the crystal above ~ 15 kcV, and thus the maximum 
cfi'cctivc area, obtained for sources at an angle of ~ 30° from the detector normal, is greater 
than the nominal 127 cm^ area of the top of the crystal at energies where the aluminum 
case is transparent. To increase the response at low energy, the aluminum case has a 3" 
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diameter beryllium window in the top surface; thus the effective area at ~ 10 keV is greatest 
for sources normal to the detector. 

The energy of a photon detected by the SDs is determined by a pulse height analyzer 
with 2752 linear channels which are then compressed to 256 quasi-logarithmic channels 
for transmission to the ground. The gain of the photomultiplier tube and the energy at 
which the lower level discriminator triggers the pulse height analyzer determine the energy 
range covered; in general a spectrum's 256 channels span two energy decades. Because 
of the priority accorded the line search, for most of the mission the gain of six of the 
eight SDs has been set to provide spectra as low as is technically feasible. As will be 
discussed below, because of detector-dependent effects, the low energy cutoff varies by 
detector. Unfortunately, after launch an electronic artifact was discovered which distorts 
the spectrum for a few channels just above the spectrum's low energy cutoff, effectively 
raising this cutoff. Although this artifact may be partially mitigated in the future by the 
calibration software, the channels of the artifact will probably never be trusted for more 
than determination of the continuum below a line candidate (Band et al. 1992). However, 
for the purpose of establishing consistency among all detectors we do search for features 
in the channels affected by this artifact when lines have been found at the same energy in 
other detectors; for example, in GB930506 the line candidate at 55 keV in one SD should 
have produced a feature just above the artifact in a second SD if it were a true absorption 
feature (Ford et al. 1994). 

After a burst triggers BATSE, the experiment collects the data necessary to analyze 
the event. In particular, the SDs accumulate a series of 192 spectra from the four most 
brightly illuminated modules. The accumulation times are varied to provide comparable 
signal-to-noise ratios (SNRs) for each spectrum, and to cover a period of 4-10 minutes over 
which the burst may last (it is of course impossible a priori to determine when the burst will 
end); during analysis on the ground various averages of these spectra are considered. These 
spectra are the basic data which are searched; later papers in this scries will consider how 
the search is conducted and how all possible averages should be treated in the statistical 
analysis. 

The mapping of the input photon spectrum into the count spectrum is traditionally 
broken into the calibration, which assigns an energy to the detector channels, and the 
detector response, which models the detector efficiency, energy broadening, and processes 
which distribute counts over a range of channels. In the BATSE SD calibration the energy 
is a monotonic function of channel number (Band et al. 1992) while the detector response 
is written as a matrix equation relating photon and channel energies. Using Monte Carlo 
simulations, the response of a BATSE SD was calculated as a function of the burst angle. 
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the angle between the detector normal and the burst (Pendleton et al. 1989); the response 
was assumed to be the same for all detectors and also independent of azimuth (the angle 
around the detector normal). The response matrix equation is 



where Cj is the count spectrum, Fj is the binned photon spectrum, and Dij is the detector 
response matrix. The energy channels need not be the same for Ci and Fj, and indeed Fj 
covers an extended energy range (particularly on the high side) to model the scattering of 
photons into the detectors' nominal energy range. The BATSE response matrix has been 
separated into direct and Earth-scattered components. The direct component includes 
photons which scatter off the spacecraft, while the Earth-scattered components consist of 
photons which impinge on the spacecraft from different directions after scattering off the 
atmosphere. The calculations included here do not include the effects of Earth-scatter. 



We use computer simulations to calculate the detectability of various types of 
absorption lines. A model photon spectrum is convolved with the appropriate detector 
response to produce a count spectrum without noise. By adding Poisson noise to this 
simulated count spectrum we create a large number of realizations which we then fit with 
continuum and continuum+line models as though they were observed spectra. Figure 1 
shows simulations of the two Ginga fines as tfiey might have appeared in the BATSE SDs. 
Line significances are calculated from the fitted parameters. We are primarily interested 
in the fraction of the simulations for which the significance of the features exceeds various 
detection thresholds; this fraction is the probability that the feature would be considered a 
detection using a given threshold. 

We map out the functional dependencies of the line detectability by varying the 
components of this procedure. First, we use response matrices corresponding to different 
burst angles. Second, we model the effect of varying the upper energy cutoff of the range 
over which the models were fit. We fix the low energy cutoff at 10 keV, a value which 
applies to some burst-detector data but is optimistic for most spectra. Third, we test the 
sensitivity to the complexity (i.e., the number of parameters) of the line model. 

Fourth, we change the signal-to-noise ratio (SNR), one of the key parameters, by 
varying the live time. The signal is the source flux convolved with the detector response. 
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3. 



METHODOLOGY 
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while the noise is assumed to result solely from the Poisson fluctuations of the signal and 

backgroTind counts (we tisc the background from an observed burst). The uncertainty 
in the background determination is not modeled and thus its contribution to the noise 
is not included in the SNR. Therefore the SNR is proportional to the square root of the 
live time. For spectra accumulated over the observed persistence times the SNR in the 
25-35 keV band is - 4.5 for GB870303 and ~ 37 for GB880205. Note that varying the SNR 
by changing the live time is not exactly equivalent to using different signal-to-background 
ratios (SBRs) since the SNR does not vary uniformly across the spectrum as the SBR is 
changed: regions of the spectrum where the signal is much larger than the background 
will see a smaller change in the SNR than regions where the background dominates. The 
SBR for om GB870303 simulation varied between 0.04 in the middle of the absorption 
line to 0.5 at ~ 150 keV, while for the GB880205 line the SBR ranged between ~ 1 over 
ii^ = 10 — 40 keV and 4.5 over E = 100 — 200 keV. However, simulations where we varied 
either the signal strength or the live time both show nearly the same dependence on SNR; 
therefore our method for calculating this dependence is robust. 

We fit the simulated spectra with the standard Marquardt-Levenberg algorithm 
minimizing (Bevington 1969, pp. 232-241; Press et al. 1992, pp. 678-683). In brief, the 
parameters of a model photon spectrum are varied to minimize the difference — quantified 
by — between the calculated count spectrum (the model photon spectrum folded through 
the instrument response) and the observed spectrum. We define with model variances, 
requiring additional terms in the gradients used by the Marquardt-Levenberg algorithm 
(see the appendix of Ford et al. 1995). 

We use the F-test to determine the hne significance. Assume the continuum-only fit 
gives Xi with ui degrees-of-freedom, while the continuum-|-line fit results in xl with 1/2. We 
define the statistic F as 



We choose a maximum probability P{> F) as the threshold for a candidate to be considered 
a line. This F-test probability is defined as the significance of the feature. Note that this 
is only the probability that the observed line is a fiuctuation; it does not consider all 
the possible "trials", the range of parameters (e.g., line centroids, widths, etc.) in which 
line-like fiuctuations could have occurred. The true probability that the line is a fiuctuation 
is the product of the F-test probability and the number of trials; therefore we require a very 
small value of the F-test probabihty to conclude the hne is real. 

We parameterize burst detectability primarily by the SNR in a number of bands across 
the spectrum; the dependence on SNR will be discussed below. Since burst continua range 
in hardness (Band et al. 1993b), the SNR in a band near the hne centroid should be used. 
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We generally characterize the BATSE spectra in the 25-35 keV band which is between the 
lines detected at 20 and 40 keV. Wc use only complete channels in these energy ranges, and 
therefore the actual energy width varies from burst to burst (the gain and consequently the 
energy-channel mapping changes with time). Since the SNR is proportional to the square 
root of the energy range AE over which it is calculated, we use 

^i^g quantity 

characterizing the signal strength in each band. 

As a consequence of statistical fluctuations, the line signiflcance in the simulated 
spectra has a broad distribution, as shown by Figure 2. The line detection probability 
p(SNR, 6) is the fraction of the simulations in which the lines were significant enough to 
be considered detections. Graphically, p(SNR, 9) is the fraction of the area under the 
curve in Figure 2 to the left of the vertical lines indicating two different possible detection 
thresholds. Since this fraction is a number based on a finite number of simulations, there 
is an associated uncertainty which we approximate as cr ~ \/p{i ~ p)/Ns, where p is the 
calculated probability and Ng is the number of simulations (usually Ng ~ 200) . For p — 
or 1 we calculate a with p = l/N^ or 1 — l/Ng, respectively. To test this approximation we 
ran 12 sets of 200 simulations at two different burst angles and two different significance 
thresholds; we found that the actual dispersion ranged from 1 to 1.4 times a. 

The line models which are used to characterize the Ginga observations are either 
additive or multiplicative Gaussians (see eqns. [3] and [4] below). Consequently, lines 
are parameterized by three quantities: an intensity (e.g., an equivalent width), a line 
centroid and an intrinsic width. At times the line models can be simplified, and fewer 
than three parameters per line are required. Models with fewer parameters may result in 
more significant fits since both the F-test and Bayesian odds compare the improvement in 
the fit to the number of added parameters. Two lines were observed from GB880205 with 
harmonically spaced energy centroids at 19.4 and 38.8 keV. Thus the spectra from this 
burst can be fit by a line model with only five parameters since the energy of the second 
line can be fixed at twice the energy of the first. The intrinsic line width is often smaller 
than, or comparable to, the instrumental resolution and the fit may not be very sensitive 
to this width. For example, the GB880205 lines have FWHM of 3.0 and 13.4 keV where 
the resolution is 5.8 and 8.8 kcV, and therefore the low energy line is unresolved while the 
high energy line will be somewhat resolved. Unresolved lines can also be modeled by a 
rectangular line profile where the fiux is zero over an energy range (the equivalent width) 
centered on the line centroid; this model requires only two parameters. For the GB880205 
lines we fit and compared the line significances in simulations with two independent three 
parameter lines (six parameters in total), two harmonically spaced three parameter lines 
(five parameters total), two independent lines with rectangular profiles (four parameters) 
and two harmonically spaced rectangular lines (three parameters). The issue is the tension 



-8- 



between the preference of statistical determinations of line significance (e.g., the F-test) for 
models with fewer parameters and the degradation in the fit when a spectrum is fit by an 
overly simplistic model. 
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4. RESULTS 



We use the spectra from the two Ginga detections to define the fine types in our 
calculations. The first is based on the SI segment of GB870303 (Graziani et al. 1992); we 
use a new fit provided by E. Fenimore (private communication, 1994), 



N{E) =5.76x10" 



E 

Too 



-1.54 



exp 



-42.5 exp 



E-2Q.Q 
2^5 



ph-cm-2-s-^-keV-\ (3) 



where E is in keV. The line was observed to persist for 4s. The second line type models the 
lines in GB880205 (Murakami et al. 1988). Again, we use a new fit provided by E. Fenimore 
(private communication, 1994), 

0.7037 r 1/^-19.4^21 
N{E) = Nc{E) - exp 
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where Nc[E) = 0.0819(^/100)-°-^^^ 
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exp 
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ph-cm 
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(4) 



V25o; 

E > 282 keV 



E < 282 keV 



note that we have expressed the Gaussians in terms of "standard deviations" a (hence the 
numerical factors). This line persisted 5.5s. The continuum we use for GB880205 does not 
correspond to the best fit Ginga models above 100 keV, but is instead the four parameter 
functional form Nc{E) (sometimes called the GRB model) which we find successfully 
describes burst spectra (Band et al. 1993b). 

For each fine type we ran simulations at burst angles of ^ = 0°, 15°, 30°, 46°, 60° and 
75° for a variety of live times which were then translated into SNR/AE'^/^; Figure 3 shows 
the results of these simulations. The spectra in GB870303 were fit over £■ = 10 — 100 keV 
(Fig. 3a) and ^ = 10 - 1200 keV (Fig. 3b), while the two lines in GB880205 were fit over 
E = \{) — 1200 kcV with a variety of line models (Figs. 3c-f). In addition, we ran simulations 
every 5° for the time the features were observed to persist (4s for GB870303 and 5.5s for 
GB880205); the resulting dependencies are shown by Figure 4. 



The results of the simulations at a given burst angle were modeled empirically by 

's - Sq 



P 



1 — tanh 



As 



/2 



(5) 



where p is the detection probability and s is the SNR/AE'^/^ (as described in §3). At s — Sq 
we have p — 1/2, while As characterizes the SNR/AE'^/^ range over which p changes from 
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to 1. Here we use an F-test probability of less than 10~^ as the detection criterion, and 
the SNR/Ai?^/^ in the 25-35 keV band. The resulting fits to the simulations are shown by 
Figure 3 and listed in Table 1. As can be seen, this model works remarkably well. It should 
be noted that these curves are fit to a small number of empirical data points, each of which 
is the average of a finite number of simulations. Consequently there is an uncertainty in 
the curves' true shape. We assume there is actually a monotonic increase in SNR/AE"^/^ 
with burst angle to achieve a given detection probability, yet the curves occasionally cross 
at small or large probabilities because of these uncertainties. 

Figure 3 and Tabic 1 establish that there is a significant dependence of the detection 
probability not only on SNR/AE'^/^ but also on burst angle. Indeed the SNR/Aii^^/^ 
necessary for a detection increases monotonically with burst angle; at small detection rates 
the empirical fits cross in a few cases, most fikely because the SNR/A£'^/^ dependence is 
undersampled, and the points were calculated with a finite number of simulations. This 
dependence on burst angle can also be seen from Figure 4 where we show the detection 
probability for the two Ginga lines with the observed line persistence times. We use 
detection criteria of F-test probabilities less than 10~^ and 10~^. 

As Figures 3c-f and the fistings for GB880205 in Table 1 also show, reducing the 
number of parameters by using a simpler line model does not change the calculated 
detectability appreciably. On the other hand, fitting the burst spectrum over a broader 
energy range makes lines more detectable; thus a line will be more significant when the 
spectrum is fit over a more extended energy range, even when the continuum model requires 
more parameters. To demonstrate this point, we fit the same 200 realizations of the model 
for GB880205 in eqn. (4) (except that in this case the break energy in the exponential 
of the GRB model Nc is 300 keV, not 250) from 10 to 100, 300, 900 and 1200 keV with 
different continuum models; the geometric mean of the F-test probabilities are presented in 
Table 2. As can be seen, the fits become more significant as the fitted spectrum is extended 
far above the high energy line at 38.8 kcV. Note that the GRB continuum model Nc is 
"correct" since this is the model used to create the spectra. Yet the simpler, less "correct" 
continuum models give more significant line fits. Adding lines to incorrect, overly simple 
continua not only fits the lines but also compensates in part for the poor continuum model, 
and therefore the improvement in the fit when the lines are added to the continuum is 
greater, as is the apparent line significance. Thus this greater line significance is misleading. 
Of course, we do not know the correct continuum shape of real burst spectra; the GRB 
model is undoubtedly also too simple. Thus our fits may also overestimate the significance 
of line candidates. 
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5. DISCUSSION 



Using a heuristic model of a spectrum with a hne we can develop some approximate 
dependencies of (Ax^), the average improvement in of the continuum+line model over 
the continuum model, on various line parameters; these dependencies explain our results, 
and provide guidance for future scalings (Ford et al. 1993). Let the observed count 
spectrum /of,s be divided into n channels of equal width Ai?, with the line falling in n^, 
channels. Assume the observed continuum is at a constant level /c, and the line is at a 
level /l (less than /c): the line profile in this simplified example is rectangular. We now 
compare a continuum model to a continuum+line model, where once again the continuum 
is assumed to be at a constant flux, and the line has a rectangular profile. On average the 
best fit continuum+line model will be the same as the observed spectrum without noise, 
while the continuum-only fit will be j^i — [(jt- — 'nL)fc + nLfLl/n. 

The noise has zero mean and non-zero variance a^. Using model variances and Poisson 
statistics we have for the ith channel — FM,i/AAtAE, where Fm^i is the model flux 
averaged over the channel, A is the effective area of the detector. At is the hve time, and 
AE is the channel width. The observed spectrum fobs is assumed to be the model spectrum 
Fm with noise, 

fobs,i = FM,i + (JiCpi (6) 

where (pi is the normalized noise distribution: = and (0?) = 1. The continuum+line 
model is "correct" and therefore 

{^l^^) = {±ih-i^lMAf)^n (7) 



err 



while for the continuum-only model 

n 



2 \ /V^ {fobs,i ~ PM,i 



,2 



= n + 



AAtAE 



E (/c - /-) + E (/^ - fmY 



(8) 



{n - riLjfc + riLjL 

Identifying the observed line width as AEn^e — tllAE, the equivalent width of the 
hne as eW — AEnneifc — fh)/ fci and the SNR for a single continuum channel as 
SNR= [UAAtAEYl'^, we find 
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Note that since SNRoc AE^/^, (SNR)^/A£' is a function of the average continuum flux 
and not the width of the energy range; indeed we parameterize the hne detectabihty by 
SNR/AE"^/^. Since our heuristic calculation used the observed count spectrum and not 
the intrinsic photon spectrum, AEune is the observed line width which results from the 
instrumental energy resolution and the intrinsic width. 

As can be seen from eqn. (2), F oc Ax^, and we find that the negative of the logarithm 
of the probability P{> F) is approximately proportional to Ax^ for a large number of 
degrees-of-freedom. Thus the absolute value of the logarithm of the line significance 
is approximately proportional to Ax^ (i.e., — logP(> F) oc Ax^)- In our case we are 
interested in the fraction of the distribution of P(> F) which meets our detection criteria, 
and therefore the relationship between detectabihty and Ax^ is more comphcated, but 
monotonic. Prom our simulations we see that the line detectabihty is dependent on the 
SNR/A£'^/^. In our simulations we fixed the hne parameters, and therefore eW and 
AEiine were not varied. However, note that the line in GB870303 is detectable at lower 
SNR/AF^/2 than the two lines of GB880205 (compare Fig. 3b to Figs. 3c-f), which can be 
understood in terms of the equivalent widths and line widths. For the photon spectrum of 
the GB870303 line the eW and AEnne are both approximately 6 keV. The 19.4 keV line of 
GB880205 has eW = 21 keV and AEu^e = 3 keV while the 38.8 keV hne has eW = 7.3 keV 
and AEiine = 13 keV. However, the analysis above deals with the observed count spectrum, 
not the photon spectrum. Thus the GB880205 line at 19.4 keV would be unresolved, and 
the 20.4 keV hne in GB870303 would be at best partially resolved, and thus both would 
be observed to have the width of the detectors' resolution of about 6 keV. Consequently 
the line in GB870303 would have a much greater equivalent width than, and comparable 
observed width to, the ~ 20 kcV line in GB880205. Similarly, the line in GB870303 would 
have a comparable equivalent width to, but smaller observed width than, the ~ 40 kcV line 
in GB880205. Thus the dependencies in eqn. (9) explain why the single line in GB870303 is 
detectable at a lower SNR/AE^/^ ^j^g^j^ ^j^g jjj^gg GB880205. 

The simulations show that extending the energy range over which the spectrum is fit 
increases the line significance, even if the continuum must be fit by a more complicated 
continuum model. This improvement in the significance results from two effects. First, 
increasing the number of degrees-of-freedom, z/2, while holding F fixed {F is the statistic 
used for the F-test — eqn. 2) decreases P{> F) (increases the line significance); for a large 
number of parameters (e.g., 6 for fitting GB880205) the factor of two increase in 1/2 resulting 
from raising the upper cutoff Ej^ax from 100 to 1200 keV can decrease P{> F) by more than 
an order of magnitude. Note that F is inversely proportional to the reduced xl = xl/^^, 
which for a good continuum-|-line fit will be approximately 1 regardless of the value of 1^2. 
Thus if Ax^ is fixed, F will most likely be fixed. Second, including more continuum in the 
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fit forces a larger difference between the continuum and continuum+line models, even when 
the added continuum is far from the energy of the candidate line feature. As eqn. (9) shows, 
increases with the ratio of the continuum counts, (n — ni)/c, to the counts in the line 
^l/l- The larger Ax^ increases F and decreases -P(> F). The first effect is in some sense 
an artifact of the statistical method while the second effect results from more fully utilizing 
the available spectral information. 

We can calculate the effect on the significance of merely changing z/2. First we find 
for each continuum type the -F1200 which gives the value of -P;y2~2i5(> -^1200) calculated 
from the simulations for z/2 ~ 215 when E^ax = 1200 keV. Using this -F1200 wc then find 
Pu2{^ -^1200) foi' the values of 1^2 applicable to the different E^ax- Table 2 shows z/2 and the 
ratio P(> F)/ Py^{> F1200) for the different combinations of continuum type and Emax- As 
can be seen, Pu2{> -^1200) decreases less rapidly than P{> F) as U2 increases. Note that the 
COMP model is poorly constrained by a continuum up to only E^^ax — 100 keV, while the 
GRB model is poorly constrained for models to E^ax — 100 and 300 keV. 

Given the above dependencies on equivalent width and SNR/A£'^/^, the dependence on 
burst angle may seem surprising. However, the energy dependence of the detector efficiency 
changes with angle, particularly at low energy where the aluminum case and beryllium 
window have different absorption properties. Between ~ 10 and ~ 20 keV the effective area 
is a maximum along the detector normal, while at higher energies, where the aluminum 
is transparent, the effective area is a maximum at a burst angle of ~ 30° (where the Nal 
crystal's projected area is greatest). 

Our simulations reveal a large distribution of fitted line parameters and significances 
for the same input photon spectrum. Figure 2 shows that even a weak line, which would be 
considered a detection in only 15% of all realizations (with an F-test threshold of 10"'^), can 
occasionally appear to be extremely significant. This distribution of line significances is not 
related to fiuctuations in the best fit line centroids, but clearly increases with the equivalent 
width of the fitted line, as demonstrated by Figure 5. The wide dispersion in the fitted fine 
centroid and equivalent width are clear from Figures 5 and 6, respectively. Consequently, 
the fitted parameters to a line detection may not be an accurate description of the true line. 

An interesting result of our simulations is that the fine detectability is nearly the same 
whether the "correct" many-parameter line model or simpler models with fewer parameters 
are used to fit the GB880205 simulated spectra. As can be seen from the definition of the 
F statistic (eqn. [2]), the significance of a line increases with a reduction in the number of 

parameters (smaller Az/) or with an improvement of the fit (larger A^;^). The width of the 
19.4 keV line is ~ 2/3 of the energy resolution while the 38.8 keV line is ~ 50% larger than 
the resolution. Thus the fits are sensitive to the line profile, and eliminating the line width 
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degrades the fit sufficiently to balance the reduction in parameters. Surprisingly, fixing the 
ratio of the line centroids to be 2 does not increase the dctectability appreciably. Note 
that reducing the number of line parameters from 6 to 5 increases F by only a factor of 
~ 1.2 (since V2 ~ 200, the difference in V2 for the two models is only ~ 0.5%) which for line 
significances of order 10~^ increases the significance (i.e., decreases P{> F)) by a factor of 
only ~ 2.7. Because of fiuctuations the best fit centroids are not at the input values, and 
therefore models with two independent line centroids result in better fits than models where 
the lines are constrained to be harmonic. The competition between simpler models and 
better fits is always present in modeling observations; with our sequence of models there is 
little preference for simplicity even when this simplicity appears to be justified on physical 
grounds. However, in other situations, such as fitting lines which are truly unresolved, 
simpler models may result in more significant fits. 

While these detectability simulations were undertaken to support the comparison of 

the BATSE and Ginga observations, they also demonstrate BATSE's detection capabilities. 
As Figure 4 shows, the SDs could have detected the lines in GB880205 up to angles of 
~ 75° while the line in GB870303 would have been detected at small angles (less than 30°) 
in about a quarter of the observations. BATSE consists of an array of detectors, and a 
burst at any point in the sky can be observed by a number of detectors at different burst 
angles. We therefore derived the average probability that the line would be detected in at 
least one detector. 

Thus we calculated the detector array's set of burst angles for 100,000 "bursts" 
distributed randomly on the sky; these sets were binned by the cosine of the burst angles. 
The probability of a detection by at least one detector was computed for each set of burst 
angles using the detection probabilities as a function of burst angle (e.g., based on the data 
in Fig. 3). For most of the mission six SDs have been set at high gain (extending their 
sensitivity to low energies), while a pair of opposite detectors have been operated at low 
gain (effectively making them useless for line searches). However, the low energy cutoff of 
one high gain detector is only ~ 25 kcV when operated at its highest gain; this detector is 
unable to detect lines at ~ 20 kcV. Only one detector of any pair of opposite detectors will 
have a positive cosine which is necessary for the burst to fall in the detector's ficld-of-view 
(unless both detectors have a cosine of zero). We therefore model the average detection 
probability using five detectors, of which four are part of two detector pairs. For comparison 
we also calculated the probabilities for two, three and four detector pairs. 

Of necessity we use a simplified model of the detector array in this calculation. The 
true energy cutoffs may differ from 10 and 1200 keV since without automatic gain control 
the gain of each detector (and therefore the energy covered) drifts with time; also, various 
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instrumental constraints cause the gain to vary. Next we assume that the data from all the 
SDs facing the burst are available. However, burst data are transmitted to Earth from only 
four detectors based on their LAD count rates, which may be affected by Earth scatter. In 
addition, the LAD and SD axes are offset by ~ 18.5°. Thus, data may be returned from one 
or two SDs with burst angles greater than 90° instead of from SDs with angles less than 90°. 
Therefore, in reality data from less than four burst-facing SDs may be returned to Earth. 
Finally, Earth blockage was not included in the calculation since in CGRO coordinates the 
Earth's position should average out over many CGRO viewing periods. 

Figure 7 shows the distribution of the cosine of the burst angle by detector rank 
assuming 4, 5, 6 and 8 detectors covered the energy range relevant for detecting lines. The 
detectors are ranked by the order of their burst angles; thus the first rank detector has the 
smallest angle, etc. As expected, operating more detectors at high gain results in a smaller 
average value of the smallest burst angle. However, the relevant quantity is the detection 
probability for the set of burst angles. For the detector array as a whole a given burst 
must be characterized by the persistence time and not by the SNR/Ai?^/^ since a variety of 
SNR/A_E^/^ values will result for a given persistence time as a function of each detector's 
burst angle. 

We calculated the average detection probability by averaging the net detection 
probability of the 100,000 bursts in each of our detector arrays, as shown by Fig. 8. This 
net probability is the probability that the line feature would be significant in at least one 
of the SDs; the probabilities for a single detector at a given angle were used to calculate 
the net probability for the entire detector array. Again as expected, the average detection 
probability increases when more detectors cover the energy range of the line features. 
The observed persistence times of the two Ginga bursts are indicated. Thus the line in 
GB870303 would be detected by BATSE with an average probability of ~ 1/8, while the 
fines in GB880205 would almost always be detected. Note that this calculation assumes 
the lines occur in spectra with the intensities reported by Ginga. The line in GB870303 is 
detectable at longer persistence times but at smaller values of SNR/Ai?^/^ than the lines 
in GB880205 because the continuum in the GB870303 spectrum is much weaker than for 
GB880205, while the line equivalent width is greater. Thus a line witfi tfic equivalent widtfi 
observed in GB870303 in a continuum as intense as GB880205 would be spectacular. 



6. SUMMARY 
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By fitting simulated BATSE spectra of the two line sets observed by Ginga we find 
that BATSE should indeed be able to detect similar lines. These siniTilations provide us 
with the detectability as a function of a spectrum's signal-to-noise ratio, a relationship 
necessary for calculations of the number of lines which would have been detected if present 
in the BATSE bursts. We note that these are calculated instrumental capabilities. 

These simulations have also provide guidance into the analysis of line candidates. 
We find that fitting spectra over a broad energy range increases the significance of line 
candidates. Decreasing the number of parameters by using physically motivated simplified 
models (e.g., assuming multiple lines are harmonically spaced or eliminating the intrinsic 
line width because it is smaller than, or comparable to, the instrumental resolution) may 
not result in the expected increase in significance, as we found in simulations with the pair 
of fines in GB880205. These lines have widths comparable to the instrument resolution, and 
therefore it is not surprising that simpfifying the profile degrades the fit. Although the fines 
were harmonically spaced in the photon spectrum, forcing the fitted lines to be harmonic 
did not significantly improve the line detectability, presumably because fiuctuations shift 
the best fit line centroids. We conclude that the degradation in the quality of the fit 
counterbalances the reduction in the number of line parameters. 

Our simulations demonstrate that BATSE is able to detect the line in GB870303 in 
about an eighth of the possible realizations and positions on the sky relative to CGRO, 
while the lines in GB880205 are almost always detectable. For most of the mission five of 
the eight spectroscopy detectors have been operated at high gain to maximize the sensitivity 
to absorption lines; increasing the gain of two more detectors (one detector cannot be 
pushed to high enough gain) will increase the average detection probabilities but not by a 
large enough factor to qualitatively change BATSE's fine detection abilities. 

The implications of these line detectabilities for issues such as consistency between the 
BATSE and Ginga observations will be explored in future papers of this series. 

We thank E. Fenimore and P. Freeman for information regarding the Ginga 
observations. We also thank the referee for pointing out the apparent improvement in 
line significance which results from the increase in the number of degrees-of-freedom as 
the fitted energy range is extended. The work of the UCSD group is supported by NASA 
contract NAS8-36081. 



-19- 



REFERENCES 

Band, D., et al. 1992, Exp. Astr., 2, 307 

Band, D., Palmer, D. Teegarden, B., Ford, L., Briggs, M., Matteson, J., Schaefer, B., 

Paciesas, W., & Pendleton, G. 1993a, in Contributed Papers of the 23rd International 
Cosmic Ray Conference, 1-105 

Band, D., ct al. 1993b, ApJ, 413, 281 

Band, D., et al. 1994, ApJ, 434, 560 

Bevington, P. R. 1969, Data Reduction and Error Analysis for the Physical Sciences 
(McGraw Hill: New York) 

Fishman, G., et al. 1989, in Proc. Gamma Ray Observatory Science Workshop (Washington: 
NASA), 2-39 

Fishman, G., et al. 1995, in preparation 

Ford, L., et al. 1993, in Compton Gamma-Ray Observatory, AIP Conference Proceedings 
280, eds. M. Friedlander, N. Gehrels and D. J. Macomb (AIP: New York), 887 

Ford, L., et al. 1994, in Gamma-Ray Bursts, Second Workshop, AIP Conference Proceedings 
307, eds. G. J. Fishman, J. J. Brainerd, and K. Hurley (AIP: New York), 261 

Ford, L., et al. 1995, ApJ, in press 

Graziani, C., Fenimore, E. E., Murakami, T., Yoshida, A., Lamb, D. Q., Wang, J. C. L., 
& Loredo, T. J. 1992, in Proc. of the Taos Workshop on Gamma-Ray Bursts, eds. 
C. Ho., R. I. Epstein, E. E. Fenimore (Cambridge U. Press: Cambridge), 407 

Hueter, G. J. 1987, PhD Thesis, UC San Diego 

Loredo, T. J. & Lamb, D. Q. 1992, in Gamma-Ray Bursts, AIP Conference Proceedings 
265, eds. W. Paciesas and G. Fishman (AIP: New York), 414 

Martin, B. R. 1971, Statistics for Physicists (Academic Press: London) 

Mazets, E. P., et al. 1981, Nature, 290, 378 

Meegan, C. A., Fishman, G. J., Wilson, R. B., Paciesas, W. S., Pendleton, G. N., 
Horack, J. M., Brock, M. N., & Kouveliotou, C. 1992, Nature, 355, 143 

Murakami, T., et al 1988, Nature, 335, 234 

Palmer, D., ct al. 1993, in Compton Gamma-Ray Observatory, AIP Conference Proceedings 
280, eds. M. Friedlander, N. Gehrels and D. J. Macomb (AIP: New York), 892 



-20- 



Palmer, D., et al. 1994a, in Gamma-Ray Bursts, Second Workshop, AIP Conference 
Proceedings 307, eds. G. J. Fishman, J. J. Brainerd, and K. Hurley (AIP: New 
York), 247 

Palmer, D., et al. 1994b, ApJ, 433, L77 

Pendleton, G. N., Paciesas, W. S., Lestrade, J. P., Fishman, G. J., Wilson, R. B., 
& Meegan, C.A. 1989, in Proc. Gamma Ray Observatory Science Workshop 
(Washington: NASA), 4-547 

Press, W. H., Teukolsky, S. A., Vetterhng, W. T., & Flannery, B. P. 1992, Numerical 
Recipes in FORTRAN, Second Edition (Cambridge: Cambridge University Press) 

Teegarden, B. et al. 1993, in Proc. of the Compton Symposium, eds. M. Friedlander, 
N. Gehrels and D. J. Macomb (AIP: New York), 860 

Wang, J. C. L., et al. 1989, Phys.Rev.Lett, 63, 1550 

Yoshida, A., Murakami, T., Nishimura, J., Kondo, L, & Fenimore, E. E. 1991, PASJ, 43, 
L69 



This preprint was prepared with the AAS IM^jX macros v3.0. 



- 21 - 



FIGURES 

Fig. 1. — Simulated count spectra of (a) the SI line in GB870303 and (b) the lines 
in GB880205 as they might have appeared in the BATSE SDs. The curve is the photon 
spectrum convolved with the detector response at a burst angle of 30°. The data points, 
rebinned into bins half a resolution element wide, are a possible realization of this count 
spectrum assuming Poisson statistics. 

Fig. 2. — Distribution of F-tcst probabilities. Based on 2400 simulations of the 
SI line of GB870303 which persisted 4s, observed at a burst angle of 15° and fit over 
= 10 — 100 keV. The two vertical lines indicate the detection thresholds of F-test 
probabilities less than 10^^ (short dashes) and 10~^ (long dashes); the fraction of the 
distribution to the left of these lines is the detection probability. 

Fig. 3. — Detection probability vs. SNR/Ai?^/^ in the 25-35 keV band for a variety 
of burst angles. From left to right, the angles are 0°, 15°, 30°, 46°, 60° and 75°. An 
F-test probability less than 10""^ was used as the detection criterion. The data points 
are the results of the simulations and the sohd curves are empirical fits to these points. 
Panels (a) and (b) are for the SI segment of GB870303 fit over F = 10 - 100 keV and 
F = 10 - 1200 keV, respectively. Panels (c)-(f) show the detectabihty of the GB880205 
line over E — 10 — 1200 keV where: (c) rectangular line profiles are used and the two fine 
centroids are required to be harmonic, E2 = 2Ei (3 line parameters); (d) rectangular line 
profiles are used and the two line centroids arc fit independently (4 parameters); (e) the 
three parameter line model is used with harmonic lines (5 parameters); and (f) the three 
parameter line model is used with the two centroids fit independently (6 parameters). 

Fig. 4. — Detection probability vs. burst angle for the two Ginga lines at the observed 
persistence times. The upper data points are for an F-test probability less than 10~^ and 
the bottom set of points for 10~^. The panels correspond to the same cases as in Fig. 3. 

Fig. 5. — Fitted equivalent width vs. F-test probability for a set of 2400 simulations of 
the SI fine of GB870303 which persisted 4s at a burst angle of 15°. The fitted equivalent 
width has been normalized to the value for the input line. 

Fig. 6. — Distribution of fitted fine centroid energies for a set of 2400 simulations of the 
SI fine of GB870303 which persisted 4s at a burst angle of 15° fitted over F = 10 - 100 keV. 
The input line centroid is indicated by the vertical dashed line. The fitted energy was 
constrained to fall between 15 and 30 keV. 

Fig. 7. — Distribution of the cosine of the burst angle by detector rank. The panels 
show the distribution for a) two detector pairs, b) 5 detectors including two pairs, c) three 
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pairs, and d) four pairs. The detectors are ranked by the size of their burst angles. For a 
detector to be included the burst must occur in its field-of-view (non-negative cosine), and 
the gain must be high enough to observe lines at ~ 20 keV. 

Fig. 8 — Average detection probabihty vs. the time a feature persisted for (a) the SI 
hue of GB870303 and (b) the lines in GB880205. From top to bottom, the curves are 
for four detector pairs (small dashes), three pairs (dash-dot-dot), 5 detectors including 
two pairs (solid), and two pairs (long dashes). The vertical lines indicate the observed 
persistence times. The burst intensities reported by Ginga are used. 



